Asymptotic silence of generic cosmological singularities 
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In this letter we investigate the nature of generic cosmological singularities using the framework 
developed by Uggla et al. We do so by studying the past asymptotic dynamics of general vacuum 
C?2 cosmologies, models that are expected to capture the singular behavior of generic cosmologies 
with no symmetries at all. In particular, our results indicate that asymptotic silence holds, i.e., 
that particle horizons along all timelines shrink to zero for generic solutions. Moreover, we provide 
evidence that spatial derivatives become dynamically insignificant along generic timelines, and that 
the evolution into the past along such timelines is governed by an asymptotic dynamical system 
which is associated with an invariant set — the silent boundary. We also identify an attracting 
subset on the silent boundary that organizes the oscillatory dynamics of generic timelines in the 
singular regime. In addition, we discuss the dynamics associated with recurring spike formation. 
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The singularity theorems of Penrose and Hawking [lj 
state that generic cosmological models contain an ini- 
tial singularity, but do not give any information on the 
nature of this singularity. Heuristic investigations of this 
issue led Belinskii, Khalatnikov and Lifshitz [2| (BKL) to 
propose that a generic cosmological initial singularity is 
spacelike, local and oscillatory. Uggla et al. [3j (UEWE) 
reformulated Einstein's field equations (EFEs) by intro- 
ducing scale-invariant variables which have the property 
that all individual terms in EFEs become asymptotically 
bounded, for generic solutions. This made it possible to 
characterize a generic cosmological initial singularity in 
terms of specific limits. The numerical study of the pic- 
ture proposed by UEWE was initiated in Ref. 0, spe- 
cializing to Gowdy vacuum spacetimes which have a non- 
oscillatory singularity. This letter presents the results of 
the first detailed study of the oscillatory asymptotic dy- 
namics of inhomogeneous cosmologies from the dynami- 
cal systems point of view introduced in UEWE. 

Here we focus on vacuum cosmologies with an Abelian 
symmetry group G 2 with two commuting spacelike 
Killing vector fields, and the spatial topology of a 3- 
torus. This is arguably the simplest class of inhomoge- 
neous models that is expected to capture the properties 
of a generic oscillatory singularity. Numerical investiga- 
tions of G 2 spacetimes supporting the BKL proposal were 
carried out by Weaver et al. in Ref. 0, . 

UEWE used an orthonormal frame formalism and fac- 
tored out the expansion of a timelike reference congru- 
ence eo by normalizing the dynamical variables with the 
isotropic Hubble expansion rate H of eo . This yielded a 
dimensionless state vector X = (E a l )(BS, where E a l are 
the Hubble-normalized components of the spatial frame 
vectors orthogonal to e ; e a — e a l di, E a l = e a l /H. 

The approach to an initial singularity will be said to be 



asymptotically silent for timelines along which E a l — > 0, 
and asymptotically silent and local for timelines along 
which E a l — > and E a l diS — > 0; in the latter case 
E a l = defines an unphysical invariant set, the silent 
boundary. [In UEWE the concept of a "silent singular- 
ity" was defined. However, the possibility of "recurring 
spike formation," discussed below, motivates the present 
distinctions and definitions.] The evolution equations for 
S on the silent boundary are identical to EFEs for spa- 
tially self-similar (SSS) and spatially homogeneous (SH) 
models, in a symmetry adapted Hubble-normalized or- 
thonormal frame . 

Motivated by the discussion in UEWE, we conjecture 
that U~^ c , the union of the bounded vacuum SH Type-I 
(Kasner) and SH Type-II subsets on the silent boundary, 
form an attracting subset that organizes the oscillatory 
dynamics of generic timelines approaching an asymptot- 
ically silent and local vacuum singularity. 

To obtain the equations for vacuum G 2 cosmologies, 
we introduce coordinates {t, x, y\, 2/2} and an orthonor- 
mal frame: e = N^ 1 d t , e.\ = ei 1 d x + e\ 2 d Vx + ei 3 cL,, 
e 2 = e 2 2 d yi and e 3 = e 3 2 d yi + e 3 3 d y . 2 , cf. Ref. @j; 
N and e a l are functions of t and x only. For compar- 
ison with previous work, e 2 is aligned with the Killing 
vector field d yi . We choose 27r-periodic coordinates x, 
yi and y 2 , yielding a spatial 3-torus topology, and a 
temporal gauge such that the area density of the G 2 - 
orbits is given by A :— (e 2 2 e3 3 ) -1 oc e _t ; this is con- 
venient since the level sets of A give a global foliation 
for maximally globally hyperbolic vacuum G 2 cosmolo- 
gies 0, and since t — ► +00 at the singularity 0- The 
G 2 symmetry implies = e 2 (/) = e^{f) for any coordi- 
nate scalar /. Thus, only iV _1 d t and ei 1 d x act nontriv- 
ially on coordinate scalars, and hence the equations of all 
spatial frame variables except ei 1 decouple; the essential 
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Hubble- normalized variables are thus E\ x := e\ x /H and 
a subset of connection components, which depend on t 
and x only. Inserting the above restrictions into the rela- 
tions in App. 5 of UEWE yields = A a = N la = A 33 = 
E12 = U2 = U3 (a = 1, 2, 3), and the spatial frame gauge 
i?i = — £23, R2 = — £31, -R3 = 0. In addition, it is con- 
venient to define: S + := ^(£22 + ^33) = — |^n> ^- := 
^(E 22 - S33), S x := ^S 23 , E 2 := 75^31, A_ := 
-^A 2 2 and A x := ^^23. The Hubble- normalized vari- 
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ables have the following physical interpretation: E+, S_, 
£ x , £2 are shear variables for eo; U = U\ describes the 
acceleration of eo; N- and A x are spatial connection 
components that determine the spatial curvature; R a 
yields the angular velocity of the spatial frame { e a }. 
The lapse function is given by A — — + 
The deceleration parameter q and the spatial Hubble 
gradient r are defined by (q + 1) := — H^ 1 8qH and 
r := — H^ 1 d\H, respectively, with do := —2(1 + £ + ) d t 
and d\ := Ex d x . These definitions yield the integrabil- 
ity condition dor — d\q = (q + 2S + ) r — (r — U) (q + 1). 

Imposing the above restrictions and gauge choices on 
EFEs in vacuum yields the following evolution equations 
and constraints: 
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doE^ = (q + 2Z + )E 1 1 
9o(l + E + ) 

do £2 
a E_ + dx A x 



(g - 2 - 3£+ + V3£_) £ 2 
(<? - 2) £_ + (r - J7) A x 
+ 2V3S 2 X - 2\/3A 2 - VS^I 
8 Q N X + =(q + 2S+) A x + (r - tf) £_ 
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where q := 2(£ 2 h + + £ 2 + £ 2 ) - | (S x - r + 
[/) [/. Since we are concerned with generic features, 
we restrict to the case £ + 7^ —1 (£+ = —1 yields 
the Minkowski spacetime). We use the gauge con- 
straint (|2a|l and the Codacci constraint l(2"c|) to solve 
for r and U and so obtain the reduced state vector 
X = ( J B 1 1 ,E + ,E2,S_,A x ,S x ,Ar_) = {Ei 1 ) © S. Note 
that the Gaufi constraint (|2b|) implies that the compo- 
nents of S are bounded. 



Because of the symmetry restrictions, E\ x is the only 
spatial frame variable in our state space; in the present 
context asymptotic silence is thus associated with E\ 1 — > 
0, while Ei 1 = is referred to as the silent boundary. 
Our numerical experiments, which employ the RNPL |l(J 
and CLAWPACK 11] packages with up to 2 16 spatial grid 
points on the x- interval (0, 2ir), indicate that asymptotic 
silence holds in the present G2 case for all timelines of 
a generic solution. Indeed, our numerical simulations in- 
dicate that max^i?! 1 ) decays exponentially. Moreover, 
they indicate that lim f ^ +oc (A := | \Ei x d x S\ | 2 ) = along 
generic timelines of a generic solution, i.e., generically the 
singularity is asymptotically silent and local, and hence 
in this case the asymptotic dynamics is governed by the 
equations on the silent boundary. 

On the silent boundary Ei 1 — 0, the integrability con- 
dition and Eq. (TcJ yield r 2 = -3/A 2 , while Eq. (|2d|) 
reduces to = (r — V3N X ) £2 ■ In contrast to Ref. p| , 
we are here concerned with the general case £2 7^ 0, 
for which r = V3N X and hence / = —1; in this case 
the equations on the silent boundary are identical to 
the Hubble-normalized equations of the exceptional SSS 
Type-_iVI models; see Wu p. 635]. 

Our numerical experiments suggest that, in addition to 
E^ -> 0, C := (U,r, A x , A_£ x ) -► holds for generic 
timelines of a generic solution when t — > +00. On the 
silent boundary Ei 1 = 0, C — yields the Kasner and 
SH Type-II subsets which are defined by = Ei 1 = 
A_ = A x = U = r, 1 = T? + + E 2 : + S 2 + S 2 , q = 2 
and = E x x = S x = A x = U = r, q = 2(1 - A 2 ), 
respectively. 

With the present gauge choices, the Kasner subset con- 
tains a subset of equilibrium points: = E\ l = S2 = 
A x = E x = A_ = U = r, 1 = + S 2 -, the Kasner 
circle, JC, which plays an essential role for the asymp- 
totic dynamics. A linear stability analysis of K. shows 



that all variables are stable when t 
(A_, E x , S2) which obey: 

7V_ = A_ e-I 1 -^)]' 
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Ei 1 and C decay exponentially and uniformly [A_£ x oc 
exp(— t)]. Here "hatted" variables are functions of x 
only. On K, S+ = t+, S_ = S_, and k(x) := 
-V3S_(o:)/[l + t+(x)]. Thus, N-, S x , S 2 are unsta- 
ble when k{x) > 1, k(x) < 0, — 1 < k(x) < 3, respec- 
tively; see Fig. Ufa). The unstable mode A_ induces 
physical curvature transitions, associated with the SH 
Type-II subset on J?/ = 0, while S x and E 2 induce 
frame transitions that lead to rotations of the spatial 
frame and multiple representations of the same solution, 
see Fig. ^b); nevertheless, for the present frame choice 
it is these gauge transitions that make repeated curva- 
ture transitions possible, and hence they have indirect 
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FIG. 1: (a) Unstable variables on /C, and (b) single transition 
sets associated with N- (dotted), E x (dash-dotted) and E2 
(solid). 



FIG. 2: Projections onto the (E+E_)-plane of a state space 
orbit along the typical timeline x = 0.3 for (a) the full G2 
system, and (b) its restriction to the silent boundary. In 
both cases the orbits approach U^ c , i.e., A — > in (a), and 
(Ei\C) -> in both (a) and (b). 



physical implications. The N-, £ x and £2 transitions 
imply that k(x) changes according to the rules k — > 2 — k, 
k — > —k and k — ► (k + 3)/(fc — 1), respectively. 

The variables and equations that describe U~ &c on 
Ei 1 = and the exceptional SH Type-VI^jyg case, as 
given by Hewitt et al. [liij . are identical. As shown in 
Sec. 5 of Ref. there exist two integrals that de- 

scribe the transition orbits. Although multiple transi- 
tions are possible, single transitions increasingly domi- 
nate. However, since frame transitions constitute gauge 
effects we will not pursue this further. What is impor- 
tant physically is that the variety U~ ac induces an in- 
finite sequence of Kasner states related by SH Type- 
II curvature transitions according to the frame invari- 
ant BKL map: u — » u — 1, if u > 2, and l/(u — 1), 
if 1 < u < 2, where u is defined frame invariantly by 
det£ Q/3 = l^jY.piY,^ = 2-27u 2 (l+u) 2 /(l + u+u 2 ) 3 . 

Numerical investigations of vacuum SH Type-VFlj/g 
and SSS Type-_iVIo models indicate that generic solu- 
tions asymptotically approach U~ ac . Our investigation 
suggests that this is also true for the evolution associ- 
ated with generic timelines of the present inhomogeneous 
vacuum G2 cosmologies, since our numerical results in- 
dicate that (Ei 1 , C, A) — > when t — > +00 for a generic 
timeline, see Fig. [21 and that thus the BKL map holds 
asymptotically for such a timeline. 

Belinskii 0] expressed concern that spatial structure, 
created by the effect of different timelines going through 
transitions at different times, could cause problems for 
the BKL scenario. Numerical experiments show that this 
is not the case for generic timelines. The reason is that 
spatial structure, created by the mechanism described 
above, develops on superhorizon scales; A — > within 
the shrinking particle horizons of a generic timeline when 
t — > +00. 

Our investigations indicate that Ei 1 — > as t — > +00 
for all timelines of a generic solution, and that (C, A) — > 
for generic timelines so that U~ ac is a local past attrac- 



tor. However, there are indications that spiky features, 
closely related to the spikes in Gowdy vacuum space- 
times, form along exceptional timelines; for such time- 
lines neither C nor A has a limit. 

Recall that for the present general G2 case, the whole 
of JC is unstable with respect to at least one of N_ , £ x , 
£2; see Eqs. (j2J. As in the Gowdy case, spike formation 
is caused by the occurrence of a zero for one of these vari- 
ables at a point (i, x(t)) when the system is close to K.. 
It follows from Eqs. ifTcl) and (|2d|) that, for a generic 
smooth solution, £2 cannot cross zero and thus pro- 
duces no spikes. Spikes in £ x are "false" (gauge) spikes, 
while spikes in N_ are "true" (physical) effects which 
yield inherently inhomogeneous dynamical features in the 
Hubble- normalized Weyl curvature scalars. Linear anal- 
ysis at K shows Ei 1 oc E\ l e _t and 



Ex l d x N- oc E^ d x N- + t A_ d x k(x) 
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The state space orbits of the spatial points outside the 
particle horizon of (t,x(t)), defined by N-(t,x(t)) = 0, 
undergo curvature transitions with iV_ = 0(1) and oppo- 
site signs on either side of x(t) when k(x) > 1; since x(t) 
does not go through such a transition this leads to the 
formation of a spike. For k(x) > 2, d\N- = E^dxN^ is 
unstable on /C, and hence grows in modulus at (t,x(t)), 
which leads to a growth in modulus of £ x . Since the par- 
ticle horizon size at (t,x(t)) is of order Ei 1 , see UEWE, 
the above implies that d x N- grows to order 1/Ei 1 and 
that diN-, and thus also A, is then 0(1) at (t,x(t)). 
It therefore follows that the dynamics fails to be local. 
Moreover, one can similarly argue that C becomes 0(1). 
Since the dynamics fails to be local at (t,x(t)), it is not 
governed by the silent boundary dynamical system. Nev- 
ertheless, our investigations indicate that the asymptotic 
dynamics is quite simple, and that it is related to that 
on the silent boundary. Numerical simulations show that 
the orbit described by S(t,x(t)) during the formation 
and smoothing out of a spike is described by the map 
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FIG. 3: Projection onto the (E+E_)-plane of a state space 
orbit undergoing a spike transition, followed by E2 and E x 
induced frame transitions and another spike transition (solid). 
The combination of N- curvature and E x frame transi- 
tions corresponding to the second spike transition is shown 
(dashed). See also Ref. [H p. 152]. 



k — > 4 — k, equivalent to a sequence of local iV_-E x -7V_ 
transitions; following Ref. |l5j, we refer to this behavior 
as a spike transition, see Fig. [3] The simplicity of this 
structure suggests that there may exist an effective dy- 
namical system governing the spike transitions, playing 
a role analogous to that of the silent boundary system. 

Numerical investigations suggest that an isolated zero 
in AU may persist as t — ► +00; if true this yields an in- 
finite sequence of recurring spike transitions. Since the 
horizon scale decays exponentially, one expects x{t) to 
converge exponentially to a point lEspike; and since the 
dynamics fails to be local during spike transitions, it fol- 
lows that C and A fail to have limits along the timeline 
x = ^spikc! we refer to such a timeline as a spike time- 
line. Timelines along which (C, A) — > as t — > +00 are 
called non-spike timelines. Since the opportunities for 
new spike formation occur at increasing time intervals 
due to the fact that JC consists of equilibrium points for 
the system, and since the horizon size decreases exponen- 
tially, we conjecture that generic timelines are non-spike 
timelines. Our analysis supports the conjecture that as 
t — > +00 the Kretschmann scalar becomes unbounded, 
also along spike timelines. 

For the Gowdy case (£2 = 0), the sequence of spike 
transitions terminates when k(x(t)) reaches the inter- 
val (0,2) on the lower part of /C; see also Ref. If 
k(x(t)) reaches the interval (0,1), the spike disappears 
completely, while in the interval (1, 2) a permanent spike 
is formed for which A — > 0. Hence, for the Gowdy case, 
A — > uniformly, and thus the singularity is asymptoti- 
cally silent and local for all timelines. 

The exponentially shrinking particle horizons of the 
spike timelines cause severe numerical difficulties. An- 
other obstacle are the subsets associated with the Taub 
points T a [cf. Fig.^a)]; in the present case T\ in particu- 
lar. For SH Type-IX models, studied by Ringstrom (it| . 



the system spends a dominant portion of its time under- 
going oscillations in the vicinity of the Taub points; this 
can be expected to hold also in the present G2 case. In 
Go cosmologies, recently studied numerically by Garfin- 
kle 0| in terms of the framework of UEWE, these is- 
sues should cause formidable problems; their resolution is 
likely to constitute a major step toward a rigorous analy- 
sis of generic spacetime singularities, and an understand- 
ing of the cosmic censorship problem. 
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